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Abstract. We present the performances of a radial velocity measurement method that we developed for A-F type stars. These 
perfomances are evaluated through an extensive set of simulations, together with actual radial velocity observations of such 
stars using the ELODIE and HARPS spectrographs. We report the case of stars constant in radial velocity, the example of a 
binary detection on HD 48097 (an A2V star, with vsin; equal to 90 kms~') and a confirmation of the existence of a 3.9 Mj up 
planet orbiting around HD 120136 (Tau Boo). The instability strip problem is also discussed. We show that with this method, it 
is in principle possible to detect planets and brown dwarfs around A-F type stars, thus allowing further study of the impact of 
stellar masses on planetary system formation over a wider range of stellar masses than is currently done. 
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1. Introduction 

Since the discovery of the first exoplanet around a solar-like 
star a decade ago ( May or & Queloz 1995| l, more than 150 plan- 
ets have been found by radial velocity surveys 1 . These sur- 
veys focus on late type stars F8) as these stars exhibit nu- 
merous lines with low rotational broadening. General charac- 
teristics of planet masses, distances to star, eccentricities (see 
e.g. |Udry et al. 2003) |Marcy et al. 2003) , as well as character- 
istics of the stars hosting those giant planets (e.g. metallicity; 
ISantos et al. 20 03 1 were derived and allow theoreticians to con- 
strain planetary system formation and evolution (e.g. planet mi- 
gration) around solar type stars. 

A general and fundamental question concerning planet for- 
mation is the impact of the mass of the central star on the for- 
mation and evolution process. We know that the disks around 
these different types of stars do not have the same properties 
at similar ages: TTauri disks of a few Myr appear to be less 
evolved than those around massive stars also of a few Myr such 
as HD 141569 or HR4796; this tends to show that these disks, 
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* Based on observations made with the ELODIE spectrograph at 
the Observatoire de Haute-Provence (CNRS, France) and with the 
HARPS spectrograph at La Silla Observatory (ESO, Chile) under pro- 
gram ID 073.C-0733. 

1 A list of discovered planets updated by Jean Schneider is available 
at http://www.obspm.fr/encycl/cat 1 .html 



as the parent stars, evolve more rapidly (Lagrange et al. 2004 1. 
The occurrence and time scale of planet formation have to be 
investigated and compared. 

Looking for planets around early type stars is a diffi- 
cult task. So far, the studies have been limited to giant stars 
( Sat o et al. 20031 lLovis et al. 2005> . These stars have small ro- 
tational velocities, but a large radius resulting in minimum pos- 
sible orbital periods of the order of 100 days or slightly less. In 
a complementary way, to focus on main sequence stars allows 
us to access smaller orbital periods and to address the question 
of evolutionary time scale. These high-mass main sequence 
stars have not been investigated so far as they exhibit fewer 
lines that are generally broadened by high rotational velocities 
(typically 100 - 200 km s -1 for A-type stars; see Fig.Q. It was 
then thought that the radial velocity method could not be ap- 
plied to those objects. Indeed, the method to process the data 
and extract the Doppler information for low-mass stars with the 
cross-correlation method, is not straightforwardly applicable to 
more massive stars IGriffin et al. 20 00 1. 

A new method for radial velocity measurements was intro- 
duced a few years ago (Chelli 2000). It consists of correlating, 
in Fourier space, each spectrum of the target star and a refer- 
ence spectrum specific to that star (built e.g. by summing all 
the available spectra of the star). This method had been applied 
to solar- type stars; we have adapted it for earlier type stars 
(Section 2). We then performed radial velocity observations 
of A-F type stars, using the fiber-fed echelle spectrographs 
ELODIE (Bar anne etal. 19961 . mounted on the 1.93-m tele- 
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Fig. 1. Examples of spectra acquired with ELODIE. Note the effect of the spectral type and vsini on the number, depth and width 
of the lines. 



scope at the Observatoire de Haute Provence (CNRS, France) 
in the northern hemisphere, and HARPS ( |Pepe et al. 20 02 1, 
recently installed on the 3.6-m ESO telescope at La Silla 
Observatory (ESO, Chile) in the southern hemisphere. 

We present here the performances obtained by applying this 
method to these spectroscopic observations of A-F main se- 
quence stars. We demonstrate in particular that it should be pos- 
sible to detect planets and brown dwarfs around A-F type stars. 
Results obtained with simulations are provided in Section 3. 
We confirm in Section 4 the accuracy of the computed ra- 
dial velocities and corresponding uncertainties in real cases, 
with examples of stars constant in radial velocity, the case of 
a binary detection and the confirmation of a 3.9 Mj up planet 
orbiting around HD 120136 (Tau Boo). Furthermore, we dis- 
cuss the first trends observed in radial velocity variations as a 
function of the spectral type in the range [A0-F7]. Finally, we 
present the uncertainties achieved for all stars already observed 
(Section 5), and the corresponding mass detection limits in the 
frame of planet and brown dwarf searches (Section 6). 



velocity associated with S (A.), we derive the cross spectrum: 
T(v) = S r (vjS*(v) = e 2 ^? |5 r (v)| 2 

The radial velocity is contained in the phase of I(v). We then 
look for the velocity V which minimizes in a least square sense 
the imaginary part of the quantity: C(v) = e~ 2 " ryA °~ I(v). The 
quantity to be minimized becomes: 
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which is minimum when V reaches the radial velocity shift U. 
Vj are discrete frequencies. The photonic radial velocity uncer- 
tainty (fRv-ph) is simultaneously calculated using the photon 
noise statistic. See Chelli (2000) for further details. 

Working in the Fourier space allows us to apply frequency 
cuts, reducing the impact of the noise (high frequencies) and of 
the variations of the continuum (low frequencies) due to stellar 
phenomena or instrumental effects. This is particularly inter- 
esting in the case of A-F type stars. 



2. Principle of the method 

To compute the radial velocity, we use the method described 
in Chelli (2000). Considering a reference spectrum S r (A) and a 
Doppler shifted one, S (A) = S , (A - A—), where U is the radial 



3. Radial velocity measurements and uncertainties 
estimates: simulations 

We first present simulations performed to test our radial ve- 
locity and photon noise uncertainty (fRv-ph) determination on 
individual measurements. 
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Fig. 2. Top: Average computed radial velocity obtained by simulations on spectra acquired with ELODIE (left) and HARPS 
(right). Error bars from the average of 100 test spectra. Bottom: radial velocity dispersions divided by uncertainties obtained 
from simulations on ELODIE (left) and HARPS (right) spectra. Conventions for symbols: vsinz < 20 kms _I (o), 20 kms -1 < 
vsinj < 70kms _1 (x), 70kms~' < vsinz' < 130 kms -1 (#), 130 kms 1 < vsin/ (+). 



3.1. Test spectra 

We consider here the case of identical spectra shifted with 
a constant radial velocity. We first consider a spectrum of a 
given star obtained with ELODIE or HARPS; this spectrum is 
smoothed then duplicated into several spectra (here, 101). We 
add noise to one of these spectra at a level corresponding to 
a typical reference spectrum, namely corresponding to the sum 
of 100 spectra with S/N equal to 200 in the case of ELODIE, and 
400 in the case of HARPS. We add noise to the other 100 spec- 
tra at a level corresponding to a typical measurement i.e. S/N = 
200 with ELODIE and 400 with HARPS. Then, these spectra are 
shifted in radial velocity with a given value (here, 300 ms _1 ), 
typically induced by the presence of a planet. 

Our method is then used to measure the radial velocities 
and the corresponding uncertainties. For the considered star, 
we obtain a distribution of computed radial velocities, charac- 
terized with its average and dispersion (if assumed gaussian). 
The advantages of these tests are: the original spectrum corre- 
sponds to a real case in terms of number and depth of spec- 
tral lines, taking into account rotational broadening and spec- 



tral type which are key parameters; the average radial velocity 
measured can be compared to the shift applied, in order to test 
the radial velocity measurement process; no effect other than 
velocity shift changes the spectra. 

3.2. Results 

The radial velocity measurements obtained are shown in Fig.|2] 
(top), for different spectral types (ranging between A0V and 
F7V), and different rotational velocities (vsin/). The measure- 
ments correspond to the average on the distribution of com- 
puted radial velocities for a given star. Error bars correspond to 
^= = e^ v _ Ph , as they correspond to the 100 test spectra for a 
given star. The average computed radial velocity is consistent 
with the initial shift RV Q equal to 300 ms _I , given the error 
bars. This demonstrates the accuracy of the computed radial 
velocities (no systematic error), in the case of identical spectra 
only shifted in radial velocity. 

The associated uncertainties are shown in Fig. [2] In all 
cases, the radial velocity dispersion <x RV on the radial velocity 
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Fig. 3. Simulated radial velocity uncertainties versus S/N (left), versus the number of spectra used to build the reference spectrum 
(center), and radial velocities obtained as a function of this number of spectra (right), in the case of an A7V star, v sin i - 90 
kms" 1 , with ELODIE. 
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Fig. 4. Measured dispersions divided by uncertainties obtained with ELODIE (left) and HARPS (right), with at least 5 measure- 
ments. S/N are those obtained during observations. Symbols are the same as in Fig.|2] 



distribution obtained for each star is consistent with the com- 
puted uncertainties eRv-ph: the distribution of J^* on the re- 
sults obtained for the different stars is centered on 0.95 ± 0.09 
for ELODIE, and on 0.94 ± 0.1 for HARPS. This demonstrates 
the quality of the photon noise uncertainty estimates. 

The same kind of tests as above (creation of 100 test spec- 
tra shifted by 300 ms" 1 , and a reference spectrum) have been 
performed to check the dependence of e^v-ph on the number n 
of spectra used to compute the reference spectrum (with S /N 
fixed to 200) and on S /N (with n fixed to 100). We find that the 
radial velocity uncertainties behave following these relations: 
ekv-Ph x where S/N is the signal-to-noise per pixel for 
one spectrum. This is useful when planning measurements, re- 
garding the exposure time necessary to reach a given eRv-p n ; 

fRV-Ph 00 -y^p" is nearly constant for n > 10. 

The results of the computation are displayed in Fig. [3] as 
dots: they are in good agreement with fits corresponding to 
these relations. Moreover, the radial velocities are not affected 
by systematic errors, even if the number of spectra used to build 
the reference spectrum is small (Fig.|5|l. 



4. Radial velocity dispersions and uncertainties: 
real case 

Our aim is to confirm the accuracy of the computed radial ve- 
locities and corresponding uncertainties in real cases. We use 
here the data available so far on our sample of stars surveyed 
with ELODIE and HARPS. By January 2005, 45 A-F type stars 
were observed at least 5 times with ELODIE (S/N equal to 196 
on average) and 13 with HARPS (S/N equal to 302 on average). 

The dispersions obtained, compared to the observed radial 
velocity uncertainties, are displayed in Fig. 0] They are con- 
sistent with an accurate computation of the radial velocity un- 
certainties in real cases, because they verify <x RV <: e RV in 

most cases, where eRv = -^ e Rv-ph + e Rv-ins • ^ n Edition to the 
photon noise, erv takes into account instrumental (in)stability 
(eteV-W 6.5 ms -1 with ELODIE, 1 ms -1 with HARPS). 

Radial velocity uncertainties observed and simulated (nor- 
malized to the same S/N) are in good agreement (see Fig. [6ji: 
in most cases, differences are smaller than 20 %. Uncertainties 
observed seem slightly larger than the ones simulated, except 
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Table 1. Radial velocity (RV) dispersion and uncertainties (measured and normalized to the same S/N, 200) for a set of stars 
observed with ELODIE and constant in radial velocity. The normalization allows us to focus the comparison between stars only 
on spectral type and v sin i. 



star 


spectral 


v.sin(i) 


number of 


dispersion on 


measured RV 


normalized RV 




type 


kms~' 


measurements 


measured RV (m s~ l ) 


uncertainties (ra s~ l ) 


uncertainties (m s~ l ) 


HD1404 


A2V 


110 


15 


203 


280 


318 


HD 102647 


A3V 


115 


17 
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222 


HD6961 
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91 
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Fig. 5. Examples of radial velocity measurements obtained in the case of stars constant in radial velocity, given the present error 
bars. 



for stars with B-V > 0.4 and v sin i < 20 km s : spectra may be 
different from one to the other in real cases, resulting in imper- 
fections in the reference spectrum, hence larger uncertainties; 
these imperfections could be negligible for late type stars, given 
the large number and depth of the lines. 

4.1. Stars constant in radial velocity 

Examples of radial velocity measurements obtained with 
ELODIE for stars that appear to be constant in radial velocity 



(given our temporal spans) are displayed in Fig. [5] and the cor- 
responding dispersions and uncertainties are given in Tabled 
The radial velocity dispersions are consistent with the uncer- 
tainties: this confirms the quality of the estimation of the later. 

4.2. HD48097: a binary system 

Among stars variable in our radial velocity measurements, we 
report here a binary detection in HD 48097 (HIP 32104), an 
A2V, vsin/ = 90km s- 1 star, with B-V = 0.063, V = 5.21, lo- 
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Fig. 6. Radial velocity uncertainties obtained in real 
cases (eRVreai) divided by simulated radial velocity 
uncertainties(eRVsimuiated), with ELODIE (left) and HARPS 
(right). Symbols are the same as in Fig. [2] 

Table 2. ELODIE orbital solution for HD48097 and Tau Boo. 



Parameter 




HD 48097 


Tau Boo 


P 


[days] 


552±17 


3.3135±0.0014 


T 


[JD-2450000] 


3005±40 


1.27±0.03 


e 




0.10+0.03 


(fixed) 


r 


[kms 1 ] 


-0.16±0.23 


-0.0650±0.004 


CO 


[deg] 


304±29 


(fixed) 


K 


[kms- 1 ] 


5.67+0.11 


0.476±0.005 


N 




15 


7 


cr(0 - C) 


[kms 1 ] 


0.19 


0.014 


a i sin i 


[AU] 


0.286 


1.510 4 


f(m) 


[M G ] 


1.0310~ 2 


3.710~ 8 


«j| 


[M ] 


2.7 


1.2 


m 2 sin i 


[M Jup ] 


430 


3.9 


a 


[AU] 


1.9 


0.046 




0.02 


-0.08 



53070 53080 53090 

JD-2400000 [days] 

Fig. 8. ELODIE radial velocity data and orbital solutions for Tau 
Boo. Top: Phased-folded velocities. Bottom: Residuals. 

thus a late K or M type dwarf (and not a white dwarf for exam- 
ple, which should have evolved faster than the primary, whose 
mass is 2.7 M ). The dispersion of the residuals is 192 ms 4 
rms, consistent with the radial velocity uncertainties (359 m s _1 
on average, 283 ms" 1 if normalized to S/N = 200); the dif- 
ference can be due to the small number of available measure- 
ments. This confirms the accuracy of the computed radial ve- 
locities in a real case. 




53000 53200 

JD-2400000 [days] 



Fig. 7. ELODIE radial velocity data and orbital solutions for 
HD48097. Top: Radial velocities. Bottom: Residuals. 



cated at 43 pc from the Sun. The orbital parameters deduced 
from a Keplerian adjustment (Fig.0 are displayed in Table|2] 
The companion is a star with a mimimum mass of 0.43 M Q , 



4.3. HD 120136: Tau Boo, a known planet 

As another example, we show the measurements obtained on 
Tau Boo (HD 120136, HR5185), anF7V star, with B-V = 0.48, 
V = 4.50, located 15 pc from the Sun. We confirm the existence 
of a planet orbiting around this star ( B utler et al. 1997i . The or- 
bital parameters (Fig. [8) deduced from a Keplerian adjustment 
are displayed in Table |2 fixing e and w to given the small 
number of measurements. They are consistent with the values 
P = 3.3128 + 0.0002 days, e = 0.02 ± 0.02, Ki = 469 ± 5 
ms" 1 previously found. Assuming a primary mass of 1.2 M , 
the minimum mass of the companion is still 3.9 Mj up with a 
semimajor axis of 0.046 AU. 

This also confirms the accuracy of the computed radial ve- 
locities, in the case of a spectral type sufficiently late to have 
been already explored. A detection of a planet orbiting an F6V 
star can be found in Galland et al. (2005b). 



4.4. Instability strip 

It appears (Fig. |4} that the late A and early F type stars ob- 
served (B-V between 0.2 and 0.4) are often highly variable 
in radial velocity. This range of B-V actually corresponds to 
the intersection of the instability strip and the main sequence, 
where we find the pulsating 5 Scuti dHandler et al. 2 002 
|Breger etaT~2 000) and y Dor stars (Mat hias etal. 20041 
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IHandler et al. 20021 . These stars can be responsible for radial 
velocity variations with an amplitude up to several kms -1 and 
periods up to several days. Additional observations are needed 
to check for each of the stars observed with B-V in the range 
0.2-0.4 and with high radial velocity variations if they are mem- 
bers of the 5 Scuti and y Dor groups. 

The studies on the associated photometric variations 
(Eye r et al. 1997| l show that the amplitude of these variations 
can reach 8 mmag with periods of 2 days in this range of spec- 
tral types. Yet, the amplitude is less than 2 mmag for the other 
A-F type stars. The latter may thus be preferential targets if 
pulsations are not desired to pollute other radial velocity varia- 
tions. 

In particular, it is probable that finding planets with the ra- 
dial velocity method will be more difficult for stars with B-V in 
the range 0.2-0.4; ways to distinguish a planetary signal from 
that of pulsations have to be investigated. 



5. Achieved uncertainties: influence of stellar 
properties 

As radial velocity uncertainties have been demonstrated to be 
accurately computed, we can discuss their values. We consider 
here the simulated radial velocity uncertainties, as they corre- 
spond to the case of identical spectra only shifted in radial ve- 
locity (keeping in mind that they are in good agreement with 
the real case, differences being less than 20 % in most cases). 

These radial velocity uncertainties are displayed in Fig. [5] 
We can see that they depend on the spectral type of the star 
(the later the spectral type, the smaller e RV , for a given range of 
v sin i), and on its rotational velocity to a higher extent. A linear 
fit of the logarithm of the radial velocity uncertainty (in m s -1 ) 
as a function of the logarithm of v sin i (in km s -1 ) gives : 

• e RV = 0.16 xv sin Z 1 54 X with ELODIE, 

• e RV = 0.032 x vsini 1 - 50 x ffe with HARPS. 

The dependance of e R y on vsin/ to the power 1.5, 
found with these fits, is consistent with the study of 
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|Bouchy et al. ( 2001) on the fundamental photon noise limit to 
radial velocity measurements, in the case of early F type main 
sequence stars. These fits apply only in the range where lines 
are resolved, and if photon noise uncertainties are large com- 
pared to the instrumental (in)stability (this corresponds typi- 
cally to vsin/ > 15 kms 1 with ELODIE). When the lines are 
not resolved (small v sin i), we reach the instrumental limit, rep- 
resented by a horizontal line in Fig.|9j in the case of ELODIE. 
Note that the dispersion around the fits results from stars with 
different spectral types. 

The achieved uncertainties on the radial velocities appear 
to be a factor of 5±1 lower with spectra obtained with HARPS 
than with ELODIE, even for early type stars which have on av- 
erage a large v sin /. A factor of 2 comes from the S/N per pixel 
corresponding to these radial velocity uncertainties (200 with 
ELODIE versus 400 with HARPS). Another factor of 2 comes 
from an increased factor of 4 in the number of pixels per spec- 
tral element, i.e. a factor of 4 in the whole flux for a given S/N 
per pixel (ELODIE and HARPS have approximatively the same 
spectral range). For the radial velocity uncertainty, these two 
effects give a factor of 2 x 2 = 4, consistent with the factor 5+ 1 
above. 

Note that we reach ELODIE and HARPS instrumental preci- 
sion for slowly rotating stars. 

6. Mass detection limits 

In the framework of searches for low mass companions, we 
now estimate what kind of such companions can be found given 
these uncertainties. 

The mass detection limits are inferred for different orbital 
periods, given radial velocity variations of + 3 6rv, assuming a 
circular orbit. They are displayed in Fig.^| The most impor- 
tant outcomes are: 

- With ELODIE, the planetary domain can be reached for 
A type main sequence stars with vsin / up to 100 kms~' 
and orbital periods less than 10 days, or with vsin / up to 
40 kms -1 and orbital periods less than 1000 days. For late 
A type stars, the accessible range is vsin / up to 80 kms~' 
and orbital periods up to 100 days. Planetary masses can be 
detected for all F type main sequence stars. 

- With HARPS, the planetary domain is accessible for all A 
and F type stars, even with large v sin i. 

For example, with ELODIE, the mass detection limit of a 10 d 
period planet around an A5 V star with v sin / equal to 60 km s 
is 4 M Jllp . Such massive close-in planets are not unexpected if 
a proto-planetary massive disk scales with the parent star mass. 
With HARPS, this detection limit decreases to 0.7 M Jup . 

7. Conclusions 

We presented in this paper the performances of a radial ve- 
locity measurement method that we developed in the case of 
A-F type stars. Radial velocities and corresponding uncertain- 
ties are shown to be accurately computed both by simulations 
and using real cases (stars constant in radial velocity, the case 



of a binary detection and the confirmation of a known planet 
orbiting Tau Boo). 

With regard to stellar properties, the achieved uncertainties 
6rv depend on the spectral type, and above all on v sin i, to the 
power 3/2: if crv is expressed in ms _1 and vsin/ in kms -1 , 
6rv behaves typically (still with a dependance on the spectral 
type) as: 

• e RV = 0.2 X vsin/ 15 X with ELODIE, if 
vsin / > 15 kms _I ; 

• e RV = 0.03 X vsin/ 15 X with HARPS. 

In particular, we have demonstrated that it should be possi- 
ble to detect extrasolar planets and brown dwarfs around such 
A-F type stars: detection limits arrive at the planetary domain 
for most of them. Given these results, we have begun a ra- 
dial velocity survey to search for these low mass companions 
around a volume-limited sample of A-F main sequence stars, 
with ELODIE and HARPS. 
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Fig. 10. Mass detection limits for A (top) and F (bottom) type stars, using ELODIE (left) or HARPS (right). 



